Ion collection by oblique surfaces of an object 
in a transversely-flowing strongly-magnetized plasma 
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The equations governing a collisionless obliquely-flowing plasma around an ion-absorbing object 
in a strong magnetic field are shown to have an exact analytic solution even for arbitrary (two- 
dimensional) object-shape, when temperature is uniform, and diffusive transport can be ignored. 
The solution has an extremely simple geometric embodiment. It shows that the ion collection flux 
density to a convex body's surface depends only upon the orientation of the surface, and provides 
the theoretical justification and calibration of oblique 'Mach-probes'. The exponential form of this 
exact solution helps explain the approximate fit of this function to previous numerical solutions. 
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Ion collection by solid objects is key to a variety 
of phenomena including the charging of dustQ and 
spacecraft 0, and the measurement of plasma parame- 
ters by electric probes In strongly magnetized sit- 
uations, where the ion Larmor radius is much smaller 
than the object, the ion collection is highly anisotropic. 
The ion motion parallel to the magnetic field can be 
approximated as governed by free-flowing compressible 
hydrodynamic equations, while perpendicular motion is 
slower, diffusive or, as in the present work, arising from 
impressed transverse flow. For convenience the object is 
here considered to be stationary and the plasma flowing, 
but this is obviously equivalent to motion of the object 
in a stationary plasma. A remarkable theoretical result 
is derived here consisting of a direct analytic solution of 
the equations governing collisionless ion collection when 
the plasma has fixed temperature, fixed transverse veloc- 
ity, negligible transverse diffusion, and arbitrary paral- 
lel velocity in the external, uniform unperturbed region. 
This solution helps resolve remaining questions concern- 
ing the calibration of so-called Mach-probes, for measur- 
ing plasma flow, and provides analytic expressions for 
the collection flux density as a function of position on a 
(two-dimensional) object of essentially arbitrary shape. 

When, as is typical, the Debye length is short com- 
pared with other lengths in the problem, the interesting 
plasma region is quasineutral. The charged thin sheath 
surrounding the object can be ignored except in so far 
as it imposes a boundary condition upon the plasma 
'presheath' solution: the 'Bohm condition' Q]. If the 
ion Larmor radius is also negligibly small, the 'magne- 
tized presheath' boundary layer where ions are acceler- 
ated across the magnetic field can also be ignored. The 
plasma transport anisotropy arising from the magnetic 
field causes the presheath perturbation to the (uniform) 
background plasma to elongate along the magnetic field 
until cross-field divergence makes up the parallel flow to 
the object [il. Full solutions to the non-linear equations 
for diffusive transverse transport have been obtained 
numer icallyd, The scaled solution, and hence the 



collected flux, depends only upon the ratio of momentum 
diffusivity to particle diffusivity. For given ratio (usually 
unity), it yields the dependence of the collected ion flux 
density on the parallel Mach- number, which amounts to 
the calibration factor for a (parallel) Mach-probe, when 
the external perpendicular velocity is ignorable. 

To measure the perpendicular plasma velocity, when 
it is not negligible, one can use oblique collection faces, 
which show sensitivity to that flow. The calibration of 
that sensitivity was calculated |7| using an extension of 
the one-dimensional approximation of @, [^, still treat- 
ing the perpendicular transport as diffusive, but account- 
ing for the boundary-condition modification that the 
transverse flow causes. This solution can be shown to be 
a Galilean transformation of the solution for zero trans- 
verse flow 3. So it is appropriate for an oblique surface of 
effectively infinite dimension in the transverse flow direc- 
tion, but finite in the direction perpendicular to both flow 
and magnetic field. Practical Mach-probes, in contrast, 
are more often multi- faceted 'Gundestrup' type[^, [lo| . 
where many short adjacent collectors are used with dif- 
ferent orientations. So it is not obvious theoretically that 
the diffusive solution should apply. And indeed the oppo- 
site limit is probably more appropriate, where diffusion 
(and viscosity) is ignored in comparison with the trans- 
verse flow. That alternative is solved here. 

The continuity and parallel momentum equations for 
the ion fluid in steady state are 



dM ^dM dn „ 
nM^— + nM— — = , 
ay ox ox 



(1) 



where M and M±_ are respectively the parallel {x-) and 
perpendicular {y-) velocities normalized to the sound 
speed, Cs, (i.e. the Mach numbers), and n is the den- 
sity. The z-coordinate is ignorable. The perpendicu- 
lar velocity is taken as a constant (positive without 
loss of generality) which is appropriate for uniform exter- 
nal flow. This perpendicular velocity assumption is the 
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only essential difference between this problem and the 
unmagnetized quasineutral fluid-plasma problem, which 
can also be analysed analytically for supersonic flow past 
a flat plate or wedge [ill. Tl^. In the outer unperturbed 
region far from the object, n = rioo and M = M^o- The 
assumption of uniform temperature is physically reason- 
able when the pressure term mostly arises from the elec- 
trons, which are assumed Boltzmann-distributed. 

These homogeneous equations can be placed into the 
following form which displays their 'characteristics': 



A/^|- + (Af + l)^ 
oy ox 



(Inn + M) = , 
(Inn- Af) = . (2) 



The solution thus requires both 

(In n + M) = const along dx = dy{M + 1) (3) 
and 

(In n - M) = const along dx = dy{M - 1) /Af^. (4) 

These will be referred to respectively as the positive and 
negative characteristics. 




FIG. 1: Construction of the solutions at different points. Po 
is in the unperturbed region. For Pi the characteristics are 
shown. Pa is in a concave region and so its positive charac- 
teristic is not tangent at P3. A value ^ > 7r/2 such as for P4 
is not problematic. 

For definiteness, we now consider plasma that is on 
the higher- a; side (to the right) of the object. Figure 
[T] will be used for illustration. For any point in the 
plasma, a positive and a negative characteristic pass 
through it. If we track the characteristics backward and 
find that they both originate in the unperturbed plasma 
at y —00, then we know that both the conditions 
In n + M = In Uoo + Moo and In n — A/ = In n^o — Moo 
(where 00 indicates values in the unperturbed region) are 
satisfied at that point. These simultaneous equations 



have only one solution: n = noo, M = Moo- The point 
is in the unperturbed region if both its characteristics 
originate there; for example: Pq. The characteristics for 
such points are straight lines with slopes M±/{M^ i !)■ 
Next we consider the possibility that just one of the 
characteristics originates not at y = —00, but on the 
surface of the object. We notice that since the positive 
characteristic dU always has more positive dz than the 
negative ([4]), the positive characteristic is always to the 
left of the negative characteristic, when tracing backward 
from a common point. Thus the positive characteristic 
will always encounter the object first; so the character- 
istic that originates on the object is the positive one. 
On that characteristic, (Inn -I- M) is constant, but not 
equal to the unperturbed value. However, because of 
the negative characteristics, the parameters do satisfy 
Inn — Af = In — Moo both at the point under con- 
sideration (e.g. Pi), and along the full length of its pos- 
itive characteristic. The only way to satisfy these two 
requirements is that, along the positive characteristic, 
M — const and n — const. If Af ~ const, then the slope 
of the characteristic, AI±/{M + 1) is constant. It is a 
straight line. 

The remaining question is, which straight-line is the 
positive characteristic? This can be answered only by 
taking into account the condition at the boundary be- 
tween the plasma and the object. The physics of the ion 
absorption, when the object is negatively charged, is ap- 
proximated by taking ions to flow into it as if they were 
flowing into a vacuum. This requires the ions to flow 
into the surface as fast as they can, consistent with exis- 
tence of an overall solution. In other words, at the object 
boundary, Af must be as negative as possible consistent 
with the overall solution. The most negative Af is equiv- 
alent to the greatest slope-angle 9 = arctan[Af j^/ (Af -I- 1)] 
(even perhaps such that > tt/2). Thus for points close 
to the boundary, the positive characteristic selected is 
the line with largest 6 that intersects the boundary (when 
tracked backwards). Such an extremum is always tangen- 
tial to the object boundary where it intersects it. Thus 
all positive characteristics that originate on the bound- 
ary do so as tangents, and for any point in the perturbed 
plasma region the positive characteristic is that straight 
line passing through the point which has greatest 9 and 
originates as a tangent on the object. Once that line 
is determined geometrically, its slope determines Af and 
hence n at all points along it. If the steepest tangent 
angle is less than 9m — M^/{Moo + 1), then the posi- 
tive characteristic does not intersect the object, but has 
slope 9 — 9„i] and the point is in unperturbed plasma. 
This specifies the entire solution for an arbitrary-shaped 
object: 

n = noo exp(Af - Afoo), Af = Af _l cot 6* - 1 . (5) 

Because the argument concerning the boundary con- 
dition is subtle, and to connect this derivation to prior 
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more specific solutions, an alternative approach giving 
the same solution is now presented. Consider first an ob- 
ject which presents to the plasma a single plane face at 
angle 0. The beginning of this face is taken as the origin. 
[We can suppose that the face which connects to it in 
an edge at the origin is parallel to the a;-axis, although 
that is not important.] It can be shown, as a straightfor- 
ward generalization of , that this system (with an ab- 
sorbing boundary) is self-similar. The solution depends 
only on the ratio x/y. Actually two self-similar forms 
satisfy the equations: either the trivial uniform solution 
AI = const, n — const] or the non-uniform solution, 
M = Mixjy =p 1, Inn = ±M -|- const (with the upper 
sign the correct choice for our configuration). The full 
solution with the appropriate boundary conditions is ob- 
tained by matching these forms by a slope-discontinuity 
along the hue x/y = x/y\m = (A^oo + i)/M±. Thus for 
this single- faced problem, for x/y < x/y\m, 

n = nooexp(M-Moo), M = M^x/y-\, (6) 

while for x/y > x/y\mi M = Moo, n = n^o- 




FIG. 2: Geometric construction of the regions of applicability 
of the solutions for a multifaceted polygonal surface. 

Now consider a multi-faceted object, as illustrated in 
figure[2l where each (planar) facet and its leading (lower) 
vertex is numbered. For the first facet (from vertices 1 
to 2), the solution © just obtained applies, in region 1. 
The region is the unperturbed, uniform plasma, and 
the line separating regions and 1 is the matching line 
x/y — {Moo + ^)/M±. Because the y- velocity is fixed and 
positive, influence cannot propagate downwards. So the 
solution on the first facet is unaffected by the presence of 
different facets at higher y. At the trailing edge of that 
facet (vertex 2) a transition to a different solution must 
occur. 



Regard vertex 2 as the new origin, relative to which the 
coordinates are {x2,y2) = {x — P2,y — 12) (where {pi,qi) 
are the coordinates of vertex i). We can then, by just the 
same approach, construct a self-similar solution to the 
governing equations ([T]) in the region adjacent to facet 
2; that solution depends only on the ratio X2/y2- This 
new solution must be matched by a slope discontinuity 
to the solution adjacent to the first facet. There is just 
one line where that matching can be made, because on 
it continuity of M requires 

A/i = M^xi/y, - 1 = M2 = M^X2/y2 - 1, (7) 

which shows xi/yi — X2/y2- The matching line is the 
projection of facet 1 beyond its upper vertex. If facet 2 
were of lesser slope than facet 1, this match line would 
not be in the plasma region; it would be inside the object. 
In that case there would in fact be no second solution, 
the solution from facet 1 would run directly into facet 2. 

Clearly the process of matching solutions can be con- 
tinued to an arbitrary number of facets. The projec- 
tions of the facets define wedge-shaped regions in which 
the solution corresponding to each facet terminates in a 
matching line. In region i, the complete solution is 

n = noo exp(M - M^o), M = M^Xi/y, ~ 1, (8) 

where {xi ,yi) = {x—pi,y — qi). If a facet lies in a concave 
region of the object, hke facets 4 and 5 of figure [H then 
it does not possess its own plasma region. Instead the 
solution(s) of the earlier region(s) apply right up to the 
respective fractions of that facet. 

There is no difficulty with the angle of the facets ex- 
ceeding 7t/2. The solution ^ applies equally well in 
regions 7 and 8 of figure [H where Xi/yi is negative. How- 
ever, if the angle of the facets exceeds tt, this solution 
does not apply because then the plasma region would lie 
to the left of the object, and the negative characteristic 
could not reach to the surface from the right-hand side of 
the object. Normally, such facets are simply part of the 
left-hand boundary of the object and are thus governed 
by the appropriate left-hand version of eq ([8]) tracked 
round the left side of the object. If (part of a) facet 
cannot be reached from either side of the object with- 
out passing through an angle of tt, then it corresponds 
to a region that has no plasma flowing to it either from 
right or left. The density there is zero. No characteristics 
originating in plasma pass through it. Finally, plasma in 
regions such that y exceeds all points on the object are 
said to be in the 'wake'. Normally the wake begins on the 
upper-most (rear-most) point of the object. At that point 
the solutions from the right and the left collide, forming 
a shock, and the simplified equations we are considering 
cannot be satisfied. For example, the temperature can- 
not be the same in the shocked region. Since the wake 
does not affect the ion collection by the object (in so far 
as equations ([1]) apply), it is not further analysed here. 
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Obviously, the continuous curved-boundary solution 
^ is the straightforward limit of the multifaceted so- 
lution ^ as the number of facets tends to infinity. The 
discussion about facet angles, etc., applies directly to tan- 
gents for the continuous case. 

The solution ^ provides an extremely simple formula 
for the ion flux to a surface not affected by concavity. 
Adding the perpendicular and parallel components, the 
total flux density is nCs{M± cos 6 — M sm9) = ncg sin 6. 
Written as flux per unit area perpendicular to the mag- 
netic field, this is 

r,! = nooC,exp[-l - (Afoc - M^cotO)]. (9) 

This form indicates, first, that for points not in a con- 
cave region of the object, the collected flux depends only 
on the angle of the surface there, and not on the shape 
of the object at smaller y. This is an important, but not 
obvious, result for establishing Mach-probe performance. 
Second, the exponential dependence upon Afoo is within 
10% of the dependence, exp[-l - 1.1{M^ - Mj_cote)], 
that fits the diffusive solutionj^, 0]; the only differ- 
ence being to replace 1.1 with 1. This agreement, 
well within most experiments' uncertainties, is an in- 
dication that the calibration is not strongly dependent 
upon the assumed cross-field transport regime. More- 
over, the present solution helps to explain why the ex- 
ponential form is such a good fit to the prior numeri- 
cal solutions, and to give confidence that it is reason- 
able to extend the form to supersonic velocities, which in 
the present context are unproblematic. Furthermore the 
present solution shows unambiguously that leading faces, 
for which 9 < 9m receive simply the unperturbed flux 
[Fii = nooCs{M±cot9 — Moo)], while trailing faces, even 
those for which Moo — Mj_ cot 9 > 1, are governed by the 
formula ([9]). The boundary conditions at the magnetic 
presheath edge arise here naturally from the analysis of 
the quasi-neutral equations. No prior assumption about 
the boundary velocity is required, although the condition 
is, of course, consistent with prior treatments. 

In a concave region of the object, e.g. at P3, the surface 



angle (local tangent) 6*8, is smaller than the characteris- 
tic's angle 9, and the distinction must be retained. This 
leads to an enhanced ion flux, equal to equation ([9]) times 
the extra factor {M± cot 9s — M± cot 9 + 1). 

For the equations used, the present solutions are ex- 
act. The formulation is, of course, approximate. Mod- 
ifications both to the boundary conditions and to the 
equations may be appropriate in some practical appli- 
cations. Nevertheless, the present analysis captures the 
dominant physics of the situation and helps to clarify 
both the operation of oblique Mach-probes and the ion 
collection by objects of essentially arbitrary shape when 
the ion Larmor radius is negligible. 
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